Beyond the Child-Langmuir Limit 
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This paper describes a new solution formulation for fully nonlinear and unsteady planar flow 
of an electron beam in a diode. Using characteristic variables - i.e., variables that follow particle 
paths - the solution is expressed through an exact analytic, but implicit, formula for any choice of 
incoming velocity vo, electric field Eq and current Jo. For steady solutions, this approach clarifies 
the origin of the maximal current J max , derived by Child and Langmuir for vq = and by Jaffe for 
Wo > 0. The implicit formulation is used to find (1) unsteady solutions having constant incoming 
flux Jo > Jmax, which leads formation of a virtual cathode, and (2) time-periodic solutions whose 
average flux exceeds the adiabatic average of Jmax- 
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Space charge limiting (SCL) current is a fundamental 
constraint on the flow of an electron beam in a diode. 
For fixed potential difference \<j)i\ and incoming velocity 
vo, the maximal sustainable current J ma x was derived by 
Child [l| and Langmuir Q for v = and by Jaffe [U 
for vq > 0. The physics origin of the SCL effect is clear: 
electromotive force from electrons in the diode limits the 
current in the diode. If incoming current is maintained 
above this maximum, then the electron density builds up 
inside the diode and a virtual cathode develops. 

The mathematical derivation of the maximal cur- 
rent in [IHIJ is based on equations for the steady, one- 
dimensional electron flow in a diode. The authors derive 
a formula relating the current and the potential jump, 
but the analysis for vq > in [3[ is quite complicated. On 
the other hand, formation of a virtual cathode is known 
to be due to cusp formation in the electron trajectories 
[H-[{|- The basic physics and technological applications 
of SCL flows and virtual cathodes are well reviewed in 
0, HI- Extensions to more general physics and geome- 
tries have been carried out, mostly using perturbation 
methods, e.g., [9] for multi-dimensional geometries. 

This paper presents a new formulation for the complete 
solution of the one-dimensional diode equations. The so- 
lution is based on characteristics (i.e., particle paths) so 
that it is an extension of [J-|6(, but it applies to both 
steady and unsteady flow and to the fully nonlinear equa- 
tions with no approximations. 

The result is an implicit solution, analogous to the im- 
plicit solution for the inviscid Burgers equation (e.g., see 
[l0|). since velocity v, density p, electric field ip, and 
spatial position x are given in terms of characteristic 
variables s and r. From this implicit formulation, it is 
staightforward to derive the maximal current that was 
first found by [H-l| and the mathematical meaning of 
the derivation is clearer in this formulation: There is a 
single derivation for vq — and vq > 0; the steps in 
the derivation are integrals along particle paths; and the 
maximal value of the flux is interpreted as a minimal 
value for the potential as a function of the incoming elec- 
tric field (for fixed flux). In addition, the implicit solu- 
tion formulation enables construction of unsteady solu- 
tions that exhibit important properties, including singu- 



larity formation corresponding to cusp formation in the 
characteristics and formation of a virtual cathode, and 
time-periodic solutions whose average flux exceeds the 
adiabatic average of J max . 

The one-dimensional continuum equations for the flux 
of electrons in a diode are 
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in which x, t, v, 0, p are the scaled position, time, velocity, 
potential and density given by 

(x, t, v, c/>, p) = (x'/L, t'/T, v/(L/T), 0'/$, p'/R) 
<S> = {m e /q e )(L/T) 2 R = e^/(q e L 2 ). 

The primed variables are unsealed, L,T are characteris- 
tic length and time values, m e , q e are electron mass and 
charge and Sq is vacuum permittivity. The boundary 
conditions at the cathode xq — and anode x\ = d are 



(4) 



so that —d>i is the potential difference across the channel 
(the negative sign is used since (f>\ < 0). 

Consider characteristic (particle path) variables in 
which x(s,t) is the position at time t = s + r for a 
particle that entered the domain at time r and is moving 
at speed v. The defining equations for s and r are 

d s x = v (5) 
x(0,t) = (6) 

t = S + T. 

Derivatives in (x,t) and in (s,r) are related by 

d s = d t + vd x (7) 
d T = d t + (d T x)d x . 
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Denote (scaled) electric field by ip = —d x <j). Since 
c00 — P, then ip(x,t) is the total mass between and 
x, plus some boundary terms, which implies 

dtil> + vd x 1> = f""(t) (8) 

for some function /"" (the four derivatives are for nota- 
tional convenience below). Combine eq. © with eq. © 
and eq. ([5]), using eq. ([Tj) , to get the following system 

= f""(s+T) 

d s v = (9) 
d s x = v. 

The general solution for this system, using eq. ([6]), is 
0(s,r) = 9(r) + f"'( S + T) 

v{s,t) = w(t)+8{t)s + /"(s + r) (10) 

X( S ,T) = w(T) S + ±6(T) S 2 +.f'( S + T)-f'(T). 

The system (|10[) provides a new general method for solv- 
ing the unsteady diode eqs. ((J)-©. 

In eq. (fTUf /, and w are related to boundary data by 

/"" = d T tPo + Jo 

= i>o-f" 
w = v - /'. 

in which Jo = Po^o and ipo = ip(x = 0) are incoming 
flux and electric field. Specification of boundary data on 
x = d requires identification of the crossing time s — T(t) 
at which characteristics (particle paths) hit x = d; i.e. 

x(T(r),T)=d. 

The density, flux and potential satisfy (using d x = (v — 
d T x)~ 1 (d s — d T ) and d T x — at x = 0) 

p = (V - drX)- 1 ^. - drW (11) 

J = v(v — d T x)~ 1 (d s — d T )ip 
Jo = (d. - 8t)1>{0,t) 

(d s -d T )4> = -{ V -drX)i>. (12) 

Eq. (fT2)) can be integrated (using 0(0, r) = 0) to get 

0(s,t) = - [ (v~d T x)ip(s',T + s- s')ds'. (13) 
Jo 

Next consider steady solutions, which cannot depend 
on r, so that /"" = Jo is a constant, and the resulting 
solutions of system ([S]) are 

ip(s) = O + Jos 

v(s) = v a + o s + JqS 2 /2 

x{s) = v s + tp s 2 /2 + J s 3 /6 (14) 

m = \{vi-v{ S f). 



In particular, the value 0i of the potential at x = d is 

0i = \{vl-v{Tf) = I^-I( Vo+ V- T+ij T 2 ) 2 . (15) 

This is equivalent to the solutions derived by Child 0- 
3] . They showed that there is maximal value Jmax of the 
current for given values of the potential difference — 0i, 
or equivalently that there is a minimal value —4> m in of the 
potential difference — 0i for given values of the incoming 
current Jo. 

The characteristic solution eq. (fl4"|) shows that there 
is a solution of the diode equations for any choice of the 
mathematically natural boundary data vq, Jo and 0o- 
The reason for the minims! potential jump — (ftmin 

(or 

equivalently the maximal current J ma x) is just that the 
function —<pi has a minimum value as 0o is varied, for 
fixed values of velocity vq and current Jo- 

In order to find this minimum value, first calculate 
d^ T by differentiating the equation vqT + 0oT 2 /2 + 
JoT 3 /6 = d and eq. (|T5j) with respect to -00 to get 

V = -^(T)- X T 2 

9^,01 = -T(v + -tp T). 

The minimal value of — 0i occurs when d^ <fii = which 
implies 

0'o = -Iv^T- 1 
d = JT 3 /6. 

At this value of T, the potential difference — 0i = — m m 
and current Jq = Jmax are 



0™„ = - W o(36d 2 J ) 1 / 3 --(36d 2 J ) 2/3 
2 o 




which is equivalent to the results of [IH1]- 

The allowable range of the parameter 0o for steady 
solutions is determined from the following two require- 
ments: The potential difference is negative (i.e., 0i < 0) 
so that electrons are driven to the anodes from the cath- 
ode; and the velocity is always positive (i.e., v(s) > for 
< s < T) since otherwise the particle paths are cross- 
ing and the model breaks down. These imply Jo > 0, 
Vq > and —JqT/2 < 0o- Stability analysis is eas- 
ily performed using the implicit solution and shows that 
the steady solution is stable for O > — 2v /T, unsta- 
ble for —2vq/T > 0o > — \/2 Jo^o and nonexistent for 
— -02 Jo^o > 0o- These results are consistent with, but 
more easily stated than those of ||. 

Next we construct several unsteady solutions that ex- 
hibit interesting behavior. As a first example, consider 
unsteady solution having constant incoming velocity Vq 
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and flux Jq. The implicit solution then has the form 

i^(s,t) = O + J Q S + f"( S + T) 

v(s,r) = v + 6 s + J ±s 2 + f"(s + T)-f"(T) 



x(s, t) 



v Q s + \e s 2 + J ±s 3 + f'(s + r) - /'(r) - /"(r) 



The corresponding potential <p is found explicitly as 

-4>(S,T) = (e +f"'(s + T))x(s,T)+ P3 ( S )-d6 

+2Jo(-/(r) + f(r + s) sf(r) \s 2 f"{r)) 

in which is defined by eq. (Tl8|) . At s = T(t), the 
equations for </>i(r) = (f>(T,r) and T(r) become 

-0i(r) - df1'+p 3 (T) + 2Jo(-f + f + -Tf'-±T 2 n 



Pl (T)-f"T + f'-f 



(17) 



in which 



pi(T) = voT + o ^ 2 + j ir 3 

j*(T) = d0o + Jo(«o^T 2 + o ir 3 + Jo^T 4 ) (18) 

/ = /(r) 
/+ = /(r + T) 

For a given function </>i, we solved the system eq. (J7| 
for /(t) and T(r) as a delay-differential equation, using 
the matlab routine ddesd, after some transformation to 
convert it into standard form for which the delays are 
backwards. This amounts to solving for incoming electric 
field ip for given values of the potential difference —(j>\. 

We present numerical computation for a solution that 
starts in the steady state F with (%, J , 4>i) = (0.5, 1, —1) 
on a system with thickness d = 4/3 (and with -0 O = 
—0.5), for t < 4. This steady state is critical in that 
the potential difference is at its minimum (i.e., — (fix = 
—<j>min) and the flux is at its maximum (i.e., Jo = Jmax)- 
The potential 4>\ varied linearly over the time interval 
4 < t < 5 up to the value (f>i = <f>i + 0.2, and then held 
constant at this value. Since this decreases the value of 
the potential jump —(f>i, it does not lead to a steady state 
solution. 

The resulting density p is presented in Figure [TJ which 
shows development of a singularity. Nevertheless, the 
function /"' remains smooth and bounded, so that im- 
plicit solution formulation remains valid up to the time of 
singularity formation. Characteristics are shown in Fig- 
ure [U which shows formation of a caustic. Note that the 
velocity becomes negative before the singularity. 

As a second example, consider unsteady solutions hav- 
ing constant incoming velocity vq but periodic flux Jo (r) 
and periodic incoming electric field iPq(t). The implicit 
form of the solution eq. ([§]) is given by 




FIG. 1: The density p for steady boundary data for which 
Jo > Jmax- Values of p are presented at 21 times starting at 
t — and at intervals of 0.4008. Near x = 0.2 value of p is 
increasing monotonically as a function of t. 




FIG. 2: Characteristics (i.e., particle paths) for steady bound- 
ary data for which Jo > Jmax- The solution breaks down 
when there is a cusp in the characteristics. 



v(s,t) = v + e(r)s + a - S 2 + f"( S + T)-f"(r) 
x(s,t) = v s + -6(t)s 2 + a o^« 3 

+f'(s+T)-f'(T)-f"(T) S 

in which ao is a constant, and 8 and / are prescribed 
periodic functions; i.e., 

/(r) - /i sm(fer/t ) 

9(t) = (9o + 6>isin(/cT/i + Ti) 

with period P = 2irto/k. The incoming flux is 

J (r) = a - e'(r). 

For given values of the constants vq, 0q, ao, ft, k, 9\ 
and n, the solution is constructed numerically: First, 
the crossing time T(t) is found by solving 



<A(s,r) = 9(t) + a s + f"'(s + t) 
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FIG. 3: The density p at various times for periodic boundary 
data for which Jo > Jmax- Values of p are presented at 21 
times starting at t = (bold curve) and at intervals of dt = 
0.1612355. The highest value of p in this figure occurs at 
t — 5dt, the time of bunching of characteristics in Figure [4] 




FIG. 4: Characteristics (i.e., particle paths) for periodic 
boundary data for which Jo > J max ■ Note that a cusp nearly 
forms in the characteristics. 



and the potential <f>i — <^(T(t),t) is found by numerical 
computation of the integral eq. (IT51) ; i.e., 



Mr) 



T(r) 



(v - d T x)i>{s',T + T(r) - s')ds' 



for a discrete set of values of r. Second, the average of 



the incoming current Jo and the adiabatic average of the 
maximal current J ma x are numerically calculated as 



Jo = P- 1 / J (t)Jt = a 
Jo 

r-P 



Jrj 



= P~ 



J maX {T){l+T\ T ))dT 



in which J ma x(''") is defined by eq. (fTlj)) using <p\ = </>i(r). 
Note that Jo is averaged over s — (i.e., x — 0) where 
dt = dr and J max is averaged over s = T(r) (i.e., x = d) 
where dt = (1 + T'(r))dr. Finally, form the difference 

Jdiff — ^0 Jmax- 

Solutions with Jdiff > (i.e., that exceed the Child- 
Langmuir limit on average) were found by Monte Carlo 
search over the values of the parameters k, 9\, t\. 
Results are shown below for [vq, ao, 0o, fi, k, 9\, t±) = 
(0.5,1,-0.5,0.1,3.8969,0.36762,2.2684) on a system 
with thickness d — 4/3. In unsealed variables, the ra- 
tio of potential energy to the energy of the entering par- 
ticles is approximately 6.5 corresponding, for example, 
to 100ey electrons on a 0.65£;V potential. The density p 
and characteristics are presented in Figures |3] and HI The 
resulting average values are Jo = 1 and J m ax — 0.85026, 
so that the average incoming current J exceeds the adia- 
batic average of the maximal current J max by an amount 
Jdiff = 0.14974, an increase of about 17%. 

Recent related work includes evidence that the average 
flux cannot exceed J m ax, under the additional constraints 
that vq — and that <pi is constant [ll[, and experimen- 
tal and numerical results for short current pulses that 
exceed J m ax, as well as for increases in current beyond 
the formation of a virtual cathode 
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The unsteady solutions constructed above suggest that 
the implicit solution formulation may be useful for ex- 
ploring additional properties of the diode equations, such 
as solutions that maximize the electric field strength and 
control methods to prevent formation of virtual cathodes. 
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